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We report results of a search for the invisible decay of the T(15) via the T(3S) — > 7r + 7r~ Y(IS') 
transition using a data sample of 2.9 fb -1 at the T(3S) resonance. The data were collected with 
the Belle detector at the KEKB asymmetric-energy e + e~ collider. No signal is found, and an upper 
limit for the branching fraction at the 90% confidence level is determined to be B(T(1S) — > invisible) 
< 2.5 x 1(T 3 . 
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An invisible particle decay mode is denned as one 
where the final state particles interact so weakly that 
they are not observable in a detector. The standard 
model (SM) predicts that there are no invisible parti- 
cles except neutrinos. If the invisible decay rate is ob- 
served to have a larger branching fraction than the SM 
prediction, it implies physics beyond the SM. One pos- 
sibility is a decay into dark matter particles, x, whose 
existence is strongly suggested by several astronomical 
observations [l[. In the SM, quarkonium decay to the 
two-neutrino final state is predicted to have a branch- 
ing fraction B(T(1S) -> vv) = (9.9 ± 0.5) x 10" 6 0. 
A much larger invisible branching fraction B(T(1S) — > 
XX) — 6 x 10 -3 is predicted [3J] for dark matter particles 
that are lighter than the b quark. Here, the pair annihila- 
tion cross section of dark matter particles to a SM quark 
pair a(xx ~ * <7<z) is estimated based on cosmological argu- 
ments, and the time-reversed reaction is assumed to have 
the same cross section, i.e. cr(qq — > xx) — f7 (xX ~~ * QQ)- 
The previous upper limits for the invisible T(IS') branch- 
ing fraction were reported by ARGUS (23 x 10~ 3 with 
90% confidence level) (3 and CLEO (50 x 10~ 3 with 
95% confidence level) [5j, which are about one order of 
magnitude above the prediction 

In this letter, we present the result of a search for the 
invisible decay of the T(IS'). The data sample used con- 
sists of 2.9fb _1 collected on the T(3S) resonance (11 



x 10 6 1(3,5)) with the Belle detector @ at the KEKB 
asymmetric-energy e + e~ (3.4 on 7.8 GeV on T(35')) col- 
lider 0]. The Belle detector is a large-solid-angle mag- 
netic spectrometer that consists of a silicon vertex de- 
tector, a central drift chamber (CDC), an array of aero- 
gel threshold Cerenkov counters (ACC), a barrel- like ar- 
rangement of timc-of- flight scintillation counters (TOF), 
and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals (ECL) located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside the coil is instrumented to detect 
K® mesons and to identify muons (KLM). 

For this search, we use the T(3S) — > 7r + 7r _ T(15) de- 
cay where only the cascade 7r + 7r~ pair is detected. If 
an T(IS') decay into an invisible final state does occur, it 
would appear as a peak at the T(IS') mass (9.46 GeV/c 2 ) 
in the distribution of recoil mass against the 7r + 7r~ sys- 
tem (M^i ) without any detected decay products from 
the Y(15). To provide a clean sample of Y(IS') decays, 
we choose to run at the T(35) resonance. Although the 
T(2S') — * 7r + 7r _ T(lS') branching fraction is about four 
times larger, the low energy of the cascade pions results 
in a low trigger efficiency. The effective cross section on 
the T(4S) resonance, a(e+e" -> T(4S) -> 7r+7r-T(lS*)), 
is much smaller than that at the T(35 l ) resonance (~ 
1/1,000 Q). To understand the reconstruction effi- 
ciency, it is essential to check that our Monte Carlo 
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(MC) simulation reproduces well the properties of the 
T(35) — > 7r + 7r~ Y(IS') transition. This is confirmed by 
using a control sample, Y(35) — * it + it~Y(1S) followed 
by T(1S) 

In our trigger system [9( , a charged track is required to 
have hits in more than half the layers of the CDC. Tracks 
that reach the outermost layer of the CDC are called 
"long tracks". To suppress beam-induced background 
events, the trigger system requires two or more charged 
tracks in an event, and at least one of them should be 
a long track. The trigger system also requires that the 
opening angle of the two tracks in the transverse plane 
be larger than a certain value for the purpose of identify- 
ing tracks being different particles. Careful evaluation of 
the trigger efficiency is important for this analysis. The 
signal events contain only two charged tracks with rela- 
tively small transverse momentum (pt)', the available en- 
ergy for the signal is 0.89 GeV, which is shared between 
the two pions. To evaluate the trigger efficiency, we took 
1.7 fb _1 of data with a trigger condition that required a 
single long track in an event with a 1/500 pre-scale rate. 

Charged tracks that are reconstructed in the CDC are 
required to originate from the interaction point: the near- 
est approach of the trajectory to the collision point is 
required to satisfy dr < 1 cm and \dz\ < 3 cm, where dz 
is measured along the direction opposite to the positron 
beam and dr in the plane perpendicular to it. We also 
require that the polar angle of the reconstructed track be 
within the detector acceptance, 17° < 6 < 150°. Charged 
kaons are distinguished from pions based on TOF, ACC 
and CDC dE/dx measurements. Electron candidates 
are identified based on the ratio of the energy detected 
in the ECL to the track momentum, the ECL shower 
shape, the energy loss in the CDC and the response of 
the ACC. Identification of muons is based on track pen- 
etration depth and the hit pattern in the KLM system. 

To reconstruct events for the control sample, we re- 
quire four charged tracks in an event, p + , fi~ , n + , and 
7r~. The T(15) — * candidates are selected by re- 

quiring 9.2 GeV/c 2 < M M+M - < 9.7 GeV/c 2 with muon 
identification for both tracks. For the pions, we reject 
tracks that are positively identified as electrons or muons. 
Figure [T] shows the distribution of the mass difference 
AM = Mu+u-Tr+Tr- — Mn+a-- The signal yield in the 
control sample is extracted from an unbinned maximum 
likelihood fit in AM. The signal shape is modeled with 
a triple Gaussian while the background shape is mod- 
eled with a first order polynomial whose slope is floated 
in the fit. We obtain 4902 ± 71 (87 ± 15) signal (back- 
ground) events in the range of 0.65 GeV/c 2 < AM < 1.15 
GeV/c 2 . Using a detection efficiency of 39.7% for the 
control sample, which is determined from MC calcula- 
tions, we estimate that 498 x 10 3 T(35) tt+tt-T (IS) 
decays are present in our data set. We compare proper- 
ties of the T(35) -> tt+tt-T(1S) decay with MC using 
events in the region of 0.875 GeV/c 2 < AM < 0.910 
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FIG. 1: Mass difference AM = M M + M - 7T+7T - -M M+M - distri- 
bution where Af„+„- lies in the T(1S) mass region. The solid 
curve shows the fits to signal plus background distribution. 
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FIG. 2: Invariant mass distribution of the two cascade pions 
for the T(3S) — > tv + tt~T(1S), T(IS') — > fi~ cascade decay 
candidates. The open histogram shows the same distribution 
for MC events. 



GeV/c 2 , where the signal purity is 99.9%. The distribu- 
tion of 7r + 7r~ invariant masses, M^+ v - is well reproduced 
by the MC as shown in Fig. [2l and the shape is consis- 
tent with CLEO results [l0( and theoretical models 11 1 . 
We also confirm that the MC reproduces well the distri- 
butions for other variables in the data such as pt of pion 
tracks, the opening angle of the pions in the transverse 
plane (tp nl7 ), the transverse momentum of the tt + tt~ sys- 
tem (p™), the polar angle of the ir + ir~ system (cos 6*^), 
and the polar angle of the muons. 

For the selection of the invisible decay candidates, 
we require two oppositely charged tracks in the event, 
i.e. 7r + 7r~ tracks. The total visible energy in the ECL 
is required to be less than 3 GeV to reject T(IS') de- 
cays into final states consisting of neutral particles. To 
minimize any possible trigger bias, we require ip„„ > 
30°, p t > 0.17 GeV/c for the tracks, and p t > 0.30 
GeV/c for at least one of the tracks. Most combi- 
natorial background is from the two-photon processes 
e + e~ — > e + e~A, where the incident e + and e~ escape de- 
tection and X — » ir + n~ , fi + p~ , n + n~ ir° or other states. 
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We reject tracks that are positively identified as electrons, 
muons or kaons. We require that no ir° candidates are 
observed in the event; we form tt° candidates using pairs 
of photons having energies greater than 20 MeV, and re- 
quiring a pair invariant mass within ±16 McV/c 2 (~ 3a) 
of the 7T° mass. Events from two-photon process are typi- 
cally boosted along the beam direction so that the vector 
Pt sum tends to be small. For further background sup- 
pression, a Fisher discriminant variable T is constructed 
with the linear combination of | cos 9^ | > Pt* an d the 
maximum energy of the 7 candidates in the event (i?™ ax ), 
T = 0.87 x I cos - 2.4 x p™ + 1.43 x £™ ax , where co- 
efficients are optimized by using MC signal events and 
background events in the Tr + n~ recoil mass sideband 
(9.42 GeV/c 2 < M^f- < 9.44 GeV/c 2 or 9.48 GeV/c 2 
< M™+ C ° U < 9.50 GeV/c 2 ). We require T < -0.7; this re- 
quirement is determined using the figure-of-merit, S/ y/B, 
where S is the number of signal events in MC and B is 
the number of background events as mentioned above. 
The overall detection efficiency, which is the product of 
the event reconstruction efficiency (9.1%) and the trig- 
ger efficiency (89.8%), is 8.2%; thus we expect 244 sig- 
nal events in our data sample for B(T(1S) — > xx) = 6 
x 10~ 3 . The event reconstruction efficiency is obtained 
from MC. 

We estimate the trigger efficiency for signal events by 
selecting two pion candidates in the single track trigger 
data under the same condition as for the signal. We ex- 
amine the trigger bits that are activated in each of the 
selected events. The track-finding efficiency of the trigger 
system is the ratio of number of signal candidates that 
satisfy the two-track trigger requirement to those that 
satisfy the single track trigger requirement. We evaluate 
the track-finding efficiency as a function of pt requiring 
a large opening angle <p„ n > 70° (the possible bias for 
two-track separation is negligible in this case). The effi- 
ciency for finding the long track is evaluated with a fur- 
ther requirement that both tracks are long tracks. The 
efficiency of the long track finding is 97.1% for p t > 0.30 
GeV/c. We also evaluate the two-track separation effi- 
ciency as a function of (p^n, where we require pt > 0.30 
GeV/c for the long track and p t > 0.17 GeV/c for the 
other track. The efficiency for the ip n7r > 30° requirement 
is 92.5%. The detection efficiency for the two-pion signal 
in the trigger system, which is the product of 97.1% and 
92.5%, is 89.8%. 

The systematic uncertainties are summarized in Ta- 
ble HI The uncertainty related to the track selection 
is estimated from the difference of the signal yield in 
data and MC for the control sample when the dr, dz, 
particle identification and other requirements are var- 
ied. The error for the 7r° veto, Fisher discriminant, 
and other requirements are also estimated in a similar 
way. Uncertainties associated with the modeling of the 
T(3S) — > n + n~T(lS) decay process are estimated by 
distorting the M w + W - distribution within the statistical 



TABLE I: Systematic uncertainties for B(T(1S) — » invisible) 
except for that from peaking background. 



Source 


{/a) 


Track selection 


. 


7T~^ veto 


2.4 


Fisher discriminant 


0.1 


Other selection requirements 


1.1 


T(35) -> 7r+7r-T(15') 


7.6 


Trigger efficiency 


8.7 


Fit bias 


0.2 


Statistics of control sample 


1.4 




2.0 


Total 


14.7 



error; possible differences in other distributions are also 
considered. The error for the trigger efficiency is obtained 
from the quadratic sum of uncertainties for the track- 
finding efficiency and the two-track separation efficiency. 
They are conservatively assigned as the differences from 
100% efficiencies. Various deviations are evaluated by 
changing the minimum p t criteria or polar angle depen- 
dences. They are negligibly small. 

For those cases where all of the T(15) decay prod- 
ucts go outside of the detector acceptance, the M^+ C °'i 
distribution still peaks at the T(15) mass and becomes 
a background to the invisible decay signal. The largest 
sources of this "peaking background" are from decays 
to oppositely charged pairs such as T(15 f ) — > or 
e + e~, where the two tracks tend to be back-to-back, so 
that when one track escapes into the forward acceptance 
hole, the other track tends to escape into the backward 
acceptance hole. The estimated contributions from dif- 
ferent decay modes, based on MC, are summarized in Ta- 
ble |n] The total number of peaking background events is 
133.2 

-lies; w here both statistical and systematic errors 
are included. 

The systematic errors for the peaking background esti- 
mation can come from uncertainties in the detection effi- 
ciency, branching fractions, MC statistics and the track- 
finding efficiency for the tracks from the T(IS') decay. 
For the contribution from the track-finding efficiency, we 
evaluate the polar angle dependence of the track-finding 
efficiency by using the control sample data. Even if we 



TABLE II: Expected number of peaking background events. 

T(1S) -> vu 0.4 ± Ol 

T(1S) -> a*V" 77.3 ± 12.0 

T(1S) -> e+e" 50.3 ± 8.2 

T(1S) -> t+t~ 5.2 ± 1.0 

Other T(1S) decay modes 0.0 + 2.8 

Other possible contributions 0.0 + 12.9 

Total 133.2 fTH 
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FIG. 3 

7r + 7r-T(lS), T(1S) -> fj, 
curve shows the fit results 



(i decay candidates. The solid 



select only one muon and two pions, we can calculate the 
momentum of the other muon track. We compare the 
data and corresponding MC. The amount of the peaking 
background differs by 3.5% between data and MC. The 
peaking background due to T(15) — > e + e~ is somewhat 
lower than that due to T(15) — > ■ This is due to the 

fact that the polar angle acceptance of the ECL (12.4° - 
155.1°) is larger than that of the CDC (17.0° - 150.0°) 
and large ECL energy deposits from an electrons can be 
present in the event, and can veto the event. Other two- 
body decays, such as T(IS') — > 7r + 7r _ and pp, are not 
observed, and their contributions are included system- 
atic error using upper limits from the PDG [13]. We do 
not observe any background contribution in MC for other 
Y(LS), T(3S I ) decays, or modes originating from initial 
state radiation. The upper limits corresponding to MC 
statistics are assigned as uncertainties. 

The signal extraction is performed by an unbinned ex- 
tended maximum likelihood fit to the M rc + coU distribu- 
tion in the range 9.40 GeV/c 2 < ikP c + c ° n < 9.52 GeV/c 2 . 
The signal shape is modeled with a double Gaussian that 
is calibrated using the control sample shown in Fig. [3] 
In the fit, the amount of peaking background is fixed 
at the estimated value and the same shape as the sig- 
nal is used. The shape of the combinatorial background 
is modeled with a first order polynomial whose slope is 
floated. Figure H shows the M^°\\ distribution. The 
extracted signal yield, 38 ± 39 events, is consistent with 
zero observed events. A x 2 test is performed for the 
M recoii distribution; we obtain x 2 /ndf = 23.4/27. The 
upper limit for the branching fraction is determined by 
the Feldman- Cousins frequentist approach [13j, taking 
into account both statistical and systematic uncertain- 
ties. We obtain B(T(1S) -> invisible) < 2.5 x 10~ 3 at 
the 90% confidence level. 

In summary, a search for the invisible decay of the 
T(15) was performed via the T(35) — > tt + tv~T(1S) tran- 
sition. In a 2.9 fb -1 data sample taken at the T(3S) 



FIG. 4: Recoil mass distribution against two pions, Af*^ 01 - . 
The solid curve shows the result of the fit to signal plus back- 
ground distributions, shaded area shows the total background 
contribution, dashed line shows the combinatorial background 
contribution, and the dot-dashed line shows the expected sig- 
nal for B(X(1S) -> XX) = 6 x 10" 3 . 



resonance, no signals were found. We obtain an up- 
per limit of 2.5 x 10~ 3 at the 90% confidence level for 
B(T(1S) — > invisible). This result disfavors the predic- 
tion in Ref. [|| for the T(15) decay to a pair of dark 
matter particles that are lighter than the b quark. 
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